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1 Introduction 

QCD sum rules [1] are very successful in determination of the masses and coupling con- 
stants of low-lying mesons and baryons. In this method a deep connection between hadron 
properties and QCD vacuum structure is established via few condensates. This approach 
is adopted and extended to many works (see for example [2] and the references therein). 
One of the important static characteristics of baryons is their magnetic moments. The 
magnetic moment of nucleon within QCD sum rules is obtained in [3, 4] using external 
field method. In [5] use is made of the QCD sum rules method in the presence of external 
electromagnetic field with field strength tensor T^ to calculate the magnetic moments of 
the baryons E c , A c containing heavy quarks. 

The goal of the present work is to calculate the magnetic moments of A b and A c in 
framework of an alternative approach to the traditional QCD sum rules, i.e., light cone 
QCD sum rules (LCQSR) (more about this method and its applications can be found in 
[6, 7] and references therein). Magnetic moments of the nucleons and decoupled baryons 
have been studied in LCQSR in [8] and [9, 10], respectively. 

The paper is organized as follows. In section 2 the LCQSR for Aq (Q = b, c) magnetic 
moment is derived. In section 3 we present numerical results. 

2 LCQSR for Aq magnetic moment 

Our starting point for determination of the Aq magnetic moment is to consider the two- 
point correlator function 

U = tJd'xe^(0\T{ VAq (x)r ]AQ (0)}\0)^ , (1) 

where T a p is the external electromagnetic field, f]\ Q is the interpolating current with Aq 
quantum numbers. It is well known that there is a continuum of choices for the baryon 
interpolating currents. The general form of Aq currents can be written as [11] 



Va q = 2 (?7 Al +tr] A , 2 ) 



(2) 



where t is an arbitrary parameter and 
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2(u T a C l5 d b )Q c + (u T a C l5 Q b )d c -(dlC l5 Q b )u c ] , (4) 



where a, b, and c are color indices. Ioffe current corresponds to the choice t — — 1. 

Let us consider phenomenological part of the correlator (2). Saturating the correlator 
with the complete set of hadron states having the same quantum numbers with Aq baryon, 
we get 



n ^ (QKIfr(Pi) ) /D , MR( „ {bMIvaJo 

n = E A _ M , (s«(pi)I^(p2)>^ vl _ M2 , (5) 



where p 2 = Pi + q, q is the photon momentum and Bi is the complete set of corresponding 
baryons having the same quantum numbers as B with masses M. 

The interpolating current couples to the baryon states with amplitudes A defined as 



°I%qI A q) = x a q u Aq (p) ■ 



(6) 



It follows from Eq. (5) that in order to calculate the phenomenological part of the correlator, 
an expression for the matrix element (B(pi)\B(p 2 )) r is needed. This matrix element can 
be parametrized as follows: 
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where f\ and f 2 are the form factors, which are functions of q 2 = (p 2 — pi) 2 and e^ is the 
polarization four vector of the photon. In the present case, in order to calculate magnetic 
moment of Aq, the values of the form factors at q 2 = are needed. Using Eqs. (5)-(7), for 
the phenomenological part of the LCQSR we get: 
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Among a number of different structures present in Eq. (8), we choose i>\4i> 2 which contains 
the magnetic moment form factor /i + f 2 . When calculated at q 2 = 0, this structure gives 
the magnetic moment of Aq baryon in units of eh/2m AQ . Isolating the phenomenological 
part of the correlator from this structure which describes the magnetic moment of the Aq 
baryon, we get 



n = -A 



A Qp 2 



ml 



"/iA c 



P\ 



ml 



(9) 



where /x Aq = (/i + /2)|ga=o- 

According to the QCD sum rules philosophy in order to construct sum rules we need 
to calculate the theoretical part of the correlator function II. Calculating correlator (1) in 
QCD we get 
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where S = CS T C, with C and T are being the charge conjugation and transpose of the 
operator, respectively. 

The perturbative contribution (i.e., photon is radiated from the freely propagating 
quarks) can easily be obtained by making the following substitution in one of the prop- 
agators in Eq. (10) 



S^^2^Jdy^y u S^ e (x-y) lfl S^ e (y)j" , (11 



where the Fock-Schwinger gauge x^A^(x) = and S^ ree is the free quark operator. In 
^-representation the propagator of the free massive quark is 

s q =^ — 7=^ '^K^i- (12) 

where tuq is the heavy quark mass and K t are the Bessel functions. Using the expansions 
for the Bessel functions 

tfi(s) ~ - + 0(x) , 

X 

K 2 (x) ~ ^2-\ + °^ 2 )^ 

and formally setting vtiq — > 0, one can obtain the well known expression of the free propa- 
gator for massless quark in x representation: 



The nonperturbative contributions can be obtained from Eq. (10) by replacing one of the 
propagators of light quarks with 
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(14) 



where A,- = |l, 75, 7 a , 2757a, a a p/\/2\ and sum over Aj is implied. 

The complete light cone expansion of the light quark propagator in external field is 
calculated in [12]. It receives contributions from the nonlocal operators qGq, qGGq, qqqq, 
where G is the gluon field strength tensor. In this work we consider operators with only one 
gluon field and neglect terms with two gluons qGGq, and four quarks qqqq, and this action 
can be justified on the basis of an expansion in conformal spin [13]. In this approximation 
heavy and massless light quark propagators read 
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(16) 



where itlq is defined from the relation 

(qWsG^a^q) = ml(qq) , 

and the operators in local part with dimension d > 5 are not taken into consideration since 
their contribution is negligible. 

It follows from Eqs. (10)-(16) that, in order to calculate the theoretical part of the 
correlator function the matrix elements of non-local operators between photon and the 
vacuum state are needed. Up to twist-4, the photon wave functions are defined in the 
following way [13]— [15]: 
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(17) 



where x is the magnetic susceptibility of the quark condensate, e q is the quark charge, the 
functions (f>(u) and ip(u) are the leading twist-2 photon wave functions, while gi(u) and 
Q2{u) are the twist-4 functions. Note that twist-3 photon wave functions are all neglected 
in further calculations since their contribution changes the results about 5%. 



The theoretical part of the correlator can be obtained by substituting photon wave 
functions and expressions of light and heavy quark propagators into Eq. (10). The sum 
rules is obtained by equating the phenomenological and theoretical parts of the correlator. 
In order to suppress the contributions of higher states and of continuum (for more details, 
see [9, 10, 16, 17]) we perform double Borel transformations of the variables p\ = p 2 and 
v\ = (p + q) 2 on both sides of the correlator. It should be mentioned here that the Borel 
transformations for K v [x) functions, which appear in the propagator of a massive quark, 
were calculated in [18]. After lengthy calculations we get the following sum rules for the 
Aq magnetic moment: 
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Here the functions ^(a,(3,x) and Fi(x) are defined as 
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where M 2 and M| are the Borel parameters. Since the initial and final baryons are the 
same, we will set M 2 = Mf = 2M 2 , from which it follows that uq = 1/2. 

It follows from Eq. (18) that the overlap amplitude Aa q needs to be known in order 
to calculate the magnetic moment. This amplitude is determined from baryon mass sum 
rules. For the mass sum rules of baryons we get (see also [19]) 
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Eqs. (20) and (21) correspond to the structures proportional to the unit operator and ]6, 
respectively. Subtraction of the continuum contribution in Eqs. (18), (20) and (21) can be 



performed by making the following substitution 

M 2 ^(a,(3,m 2 Q /M 2 ) -> ^^r- / 2 dae-/ M / *( s - tmlf'H^^it - l)- 1 ,(22) 

for a > and n > 0. 

3 Numerical analysis 

This section is devoted to the numerical analysis of the sum rules for the magnetic moment 
of Aq baryon. It follows from Eq. (18) that the main input parameter of the LCQSR 
is the photon wave function. It was shown in [13, 14] that the leading twist-2 photon 
wave functions receive minor corrections from the higher conformal spin, so that to a good 
approximation they can be described by their asymptotic forms. In further numerical 
calculations we will use the following forms of the photon wave functions [13, 15]: 

4>(u) = 6u(l — u) , i>(u) = 1 , 
gi(u) = --(1 - w)(3 - u) , g 2 (u) = --(1 - uf . 

The values of the other input parameters which we need in the numerical calculations are: 
/ = 0.028 GeV 2 , x = -4.4 GeV~ 2 [20] (in [21] this quantity is estimated to be x = 
-3.3 GeV- 2 ), (qq)(l GeV) = -(0.243) 3 GeV 3 , m 2 = (0.8 ±0.2) GeV 2 [22], m c = 1.3 GeV 
and mj = 4.8 GeV. 

In view of the fact that the magnetic moment is a physical quantity it must be inde- 
pendent of the auxiliary parameters t, Borel mass M 2 , as well as the continuum threshold 
So- Therefore calculation of the magnetic moment requires simply finding a region where 
H\ Q is practically independent of the above-mentioned parameters. 

Our strategy in trying to resolve this problem consists of three steps. In the first step, 
in order to find the working region, where /i\ Q is supposed to be independent of the Borel 
parameter M 2 , we study the dependence of the magnetic moment /j,\ q on M 2 at several 
predetermined values of the threshold sq and at three different values of t. Along these lines, 
in Fig. (1) we present the dependence of /i\ c on M 2 at so = 10 GeV 2 and so = 15 GeV 2 , 
for different values of t. Similarly, Fig. (2) depicts the dependence of fi\ b on M 2 at two 
different values of the threshold, Sq = 40 GeV 2 and s — 45 GeV 2 . From Fig. (1) we 
observe that the working region of M 2 which is consistent with our requirements lies in the 
interval 3 GeV 2 < M 2 < 5 GeV 2 for the A c case and 15 GeV 2 < M 2 < 25 GeV 2 for the 
At case. Obviously, from both these figures we observe that fX\ Q is reasonably insensitive 
to the different choices of the continuum threshold, while it seems to be sensitive to the 
parameter t. 

The next step, of course, is to explore the physical region for the parameter t. For this 
purpose we have used the mass sum rules given in Eqs. (20) and (21), both of which are 
required to be positive. Following this line of reasoning, we present in Figs. (3) and (4) 
the dependence of the mass sum rule (20) on cos#, where 9 is determined from the relation 
tan# = t, for A c and A&, respectively. These figures depict that the relevant physical 
regions for the parameter t are given by —0.78 < cos# < 0.7 for A c and —0.8 < cos# < 0.7 



for A^ respectively. The analysis of (21), clearly, leads to the conclusion that this sum 
rule is positive for arbitrary values of the parameter t. As a result of these observations, 
the physical region of t which guarantees the positiveness of the sum rules (20) and (21) 
separately, is confined to the interval —0.78 < cos9 < 0.7. Having this restriction on t, our 
final attempt is to determine the value of the magnetic moment fi\ Q . For this purpose, we 
must find a region of t where /j,\ q is independent of this parameter. 

In Figs. (5), we present the dependence of the magnetic moment of heavy baryon 
A c on cos# at the fixed value of the Borel parameter M 2 = 4 GeV 2 for two different 
choices of the continuum threshold sq = 3 GeV 2 and sq = 4 GeV 2 . Similarly, depicted in 
Fig. (6) is the dependence of the other heavy baryon A& on cos 6 at M 2 = 20 GeV 2 and 
So = 40 GeV 2 ; sq = 45 GeV 2 . We observe from these figures that the magnetic moment is 
quite stable in the region —0.25 < cos6 < +0.5, and practically seems to be independent 
of cos#, t and the continuum threshold Sq. As a result of all these considerations we obtain 
for the magnetic moment 

/i Ac = (0.40 ± 0.05) jujv , 

/xa 6 = (-0.18 ±0.05) /Mr, (23) 

where /in is the nucleon magneton. 

Finally we present a comparison of our result on fi\ Q with the existing theoretical 
calculations in literature. For the magnetic moment fi\ c the traditional QCD sum rules 
predicts fi\ c = (0.15 ± 0.05) fiN [5]. Our result on /j,\ c is close to the non-relativistic quark 
model prediction [23, 24], but there is substantial difference with the results predicted in 
[5]. In our opinion, this discrepancy can be attributed to the fact that the results presented 
in [5] were calculated for the choice t = — 1, which is unphysical in our case. Magnetic 
moments of triplet and sextet heavy baryons have been calculated using heavy hadron 
chiral perturbation theory (HHChPT) in [25, 26, 27], and A c and A& baryons that we are 
interested belong to the triplet. It was shown in these works that in HHChPT the leading 
term in magnetic moments is proportional to eQ/m c , where Q is the heavy quark. The 
corrections to the leading term appear at the order of 0(1/i71qA 2 ), where A x ~ 1 GeV, 
and four arbitrary constants are required. If we set all four arbitrary constants to zero, for 
the magnetic moments of A c and A& baryons HHChPT predicts 

fi Ac ~ 0.47 fj, N , 

/i Ab ~ -0.23 fj, N . (24) 

From Eqs. (23) and (24) we conclude that both approaches lead to close values for the 
magnetic moments of the heavy A c and A& baryons. The small difference between the 
predictions on the magnetic moments of heavy baryons of the two approaches is due to the 
subleading terms which are neglected in deriving Eq. (24). 

A measurement of the magnetic moments of heavy baryons represents an experimental 
challenge. Few groups are contemplating the possibility of performing magnetic moments 
in the near future (BTeV and SELEX) [28]. 
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Figure captions 

Fig. (1) The dependence of the magnetic moment fi\ c on M 2 at two different values of 
the continuum threshold sq = 1 GeV 2 and So = 15 GeV 2 , for several fixed values of the 
parameter t. Here in this figure and in all following figures the magnetic moments of A c 
and A?, baryons are given in units of the nucleon magneton /xjv- 

Fig. (2) The dependence of the magnetic moment fi\ b on M 2 at two different values 
of the continuum threshold sq = 40 GeV 2 and sq = 45 GeV 2 , for several fixed values of the 
parameter t. 

Fig. (3) The dependence of the mass sum rule m^ c on cos 9, at two different values 
of the continuum threshold s = 10 GeV 2 and s = 15 GeV 2 . 

Fig. (4) The dependence of the mass sum rule m\ b on cos 6, at two different values 
of the continuum threshold sq = 40 GeV 2 and so = 45 GeV 2 . 

Fig. (5) The dependence of the magnetic moment fi^ c on cos#, at M 2 = 4 GeV 2 and 
at two different values of the continuum threshold Sq = 3 GeV 2 and s — 4 GeV 2 . 

Fig. (6) The dependence of the magnetic moment fi\ b on cos#, at M 2 = 20 GeV 2 and at 
two different values of the continuum threshold Sq = 40 GeV 2 and Sq = 45 GeV 2 . 
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